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Abstract: We estimate the ratio of the thermalization time over the freeze-out time 
using a holographic AdS-Vaidya model for quark-gluon plasma formed in the heavy 
ion collisions. In the model the process of thermalization is described as formation of 
the black hole in AdS space while dethermalization process, related with the freeze- 
out, as the black hole evaporation due to the Hawking radiation that is modeled by 
the Vaidya metric with a negative mass. In this model the thermalization takes place 
only at small scales and absent in the infrared region. At small scales the system 
tends towards a state of the thermal equilibrium only for the short time after which 
the processes of dethermalization starts. In this simple model the dethermalization 
time has a low bound about 7 fm/c which is consistent with experimental data. 
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1. Introduction 



In the process of collision of two ions one observes a rapid local thermalization of 
the system of partons with further expansion and cooling of system which leads to 
hadronization and multiple production of particles. Experiments indicate to a very 
short thermalization time, Ttherm ~ 1 fm/c for the quark-gluon plasma(QGP) formed 
in heavy ion collisions, while the freeze-out time is of order 20 fm/c |L|, 0, ||, f|, [5], §]. 

In more details, the time schedule of the process is the following. Up to a time 
~ 0.02 fm/c "hard" processes take place and they are responsible for "hard" 
particles, which can be observed at detectors. Up to a time ~ 0.2 fm/c "semi-hard" 
processes take place and they produce the most of the "multiplicity" in the final 
state. Then at the thermalization time of order r^ erm ~ 1 fm/c the system reaches 
a local thermal equilibrium state, called QGP. After that the evolution of QGP is 
described by equations of hydrodynamics and after the time of order r hadr ~ 10 
fm/c, when due to the separation of the colliding ions the temperature becomes 
lower than the deconfinement temperature, a hot hadron gas is formed. Upon the 
further expansion and cooling, around the freeze-out time, we refer to this time as 
the dethermalization time Tdet ~ 20 fm/c, the density of the hadron gas become 
sufficiently low and the system decays into free hadrons, which can be observed at 
detectors. Therefore in experiments on heavy ion collisions there is the following 
hierarchy of time scales: 

Ttherm *^ Tfiydro *^ Tfiadr ^ 7~det- (I'-O 

There are many attempts to describe the process of heavy ion collisions and QGP 
formation in the QCD framework. One of difficulties is that one has to compute the 
time depended correlation functions in the strong coupling regime since the QGP 
is a strong coupling system ||. In the recent years a powerful approach to these 
problems is pursued which is based on a holographic duality between the strong 
coupling quantum field in (i-dimensional Minkowski space and classical gravity in 
d + 1-dimensional anti-de Sitter space (AdS) |7|, ||, |9|. In particular, there is a 
considerable progress in the holographic description of equilibrium QGP |10 . 



Holographic thermalization describes the thermalization process as a process of 
formation of a black hole in AdS. There are several scenarios to produce a black hole 
in AdS 0-0. One of them uses the AdS-Vaidya metric [EJ, [27], 



In this paper we consider an influence of the Hawking radiation on the ther- 
malization process. We describe the evaporation of black hole due to the Hawking 
radiation by the Vaidya metric with negative mass. We use a special form of the 



AdS-Vaidya metric, see eq. (|4.1| ) below, to describe thermalization and subsequent 
dethermalization. We show that the evaporation of the black hole in AdS leads to an 
interesting phenomena in the d- dimensional Minkowski space-time - thermalization 
is possible only at small distances and impossible in the infrared region. 
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The Hawking radiation from a Vaidya black hole was considered in [34] where it 
was shown that allowance for the nonstationarity reduces to allowing for the depen- 
dence of the temperature of the Hawking radiation on the retarded time. Note that 



the Hawking radiation in holographic approach was also discussed in [[35 

The paper is organized as follows. In Sect. 2 we remind the basic formulas 
related with the AdS- Vaidya metric which are used for the holographic thermaliza- 
Sect.3 discusses the AdS- Vaidya metric with negative mass to model 



tion 26 - 29 



the Hawking radiation. We show that thermalization for the two-point correlation 
functions holds only at small distances. In Sect. 4 we discuss attempts to fit the 
time hierarchy (LI) obtained in Pb+Pb collisions at LHC and in Au+Au collisions 
at RIHC, by a special form of the AdS- Vaidya metric. In Appendix we present 
an explicit formula for a non-equal time two points correlation function for the de- 
thermalization process. This formula is a counterpart of non-equal time two points 
correlation functions for the thermalization process found in |[26|-[|29| . 



2. Holographic Thermalization with AdS-Vadiya metric 

Holographic duality prescribes that the vacuum correlation functions in quantum 
field theory in the strong coupling regime in 4-dimensional Minkowski space can be 
computed by using the action functional for corresponding fields in d + 1- dimensional 
AdS 0, @, [J. Thermal states also admit a holographic description by using a black 
hole AdS metric || 38 1. The holographic description of equilibrium QGP was tested 



in numerical works, see the review paper |IJ and refs therein. 



Hypothesis on holographic thermalization asserts that the thermalization process 
can be described in the dual framework as the process of formation of a black hole 
in AdS. The process of formation of a black hole can be initiated by a perturbation 
of the initial AdS metric. Formations of black holes in gravitational collapse or in 
the collision of ultra-relativistic particles are traditional difficult problems in general 
relativity, see for example |37| and refs therein. In AdS there are specific features re- 



lated with the compression of AdS metric. In the context of holographic description 
of heavy-ions collisions formation of black holes in collision of gravitational waves in 
AdS has been considered in numerous works |11|, |L|, [13], [15], [I7|, [L6|, |1^, |I9~| . Gravi- 



tational collapse occurring as a result of the weak perturbation in AdS is considered 
in [^U| , where it is shown that the system evolves to the state of thermal equilibrium 
almost instantly. This result stimulates application of the AdS- Vaidya metric for the 



thin shell in the holographic description of thermalization |32], [J3[ [26], £7], |29| . 

The (d + l)-dimensional infalling matter shell in AdS in Poincare coordinates is 
described by the Vaidya metric 

ds 2 = — [-(!- m(v)z d ) dv 2 - 2dzdv + dx 2 ] , (2.1) 
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where v is the null coordinate, x = (x\, . . . , Xd-i) are the spatial coordinates on the 
boundary z = and we have set the AdS radius equal to 1. We take m(v) in the 
form 

m(v) = M6{v), (2.2) 

where M is a constant and 8(v) is the Heaviside function. 
For m(v) = M the change of variables 



dv = dt 



dz 



Mz a 



(2.3) 



brings ( |2.1| ) to the standard metric of the black hole in AdS in the Poincare coordi- 
nates 



ds 2 



1 



dz 2 



- (1 - Mz d ) dt 2 + + rfx : 

v ; 1 - Mz d 



(2.4) 



For v < the metric ( |2.1| ) with m(v) in the form ( |2.2j ) is just the AdS metric. 

The holographic prescription permits to calculate an equal-time two-point func- 
tion of operators 0(x, t) with large conformal dimension A in dual theory of gravity 
by using the geodesic approximation [ |36|j : 



(O(x,t)O(0,t)) ~exp[-A£(x,t)]. 



(2.5) 



Here £(x, t) is the length of the geodesic that begins and ends at the boundary points 
(0,£) and (x,t). The vacuum correlation functions correspond to the computation of 
geodesies in AdS, thermal ones in AdS with a black hole (BHAdS) while correlation 
functions describing the process of thermalization correspond to geodesies in the AdS- 
Vadiya ( |2.ip with (|2.2|). These geodesies are studied in papers |£| 



These considerations give an estimation of the thermalization time for two equal time 
points at the boundary. In particular, for d = 2 the thermalization time is equal to 
the half of the space distance between these points. For d > 2 the thermalization time 
decreases. The Vaidya Reissner-Nordstrom AdS metric increases the thermalization 
time in all dimensions |]30l |3T 



3. The Hawking Radiation and Dethermalization 

In this section we consider the Hawking radiation in the holographic approach by 
using the AdS- Vaidya metric and show that the Hawking radiation prevents to ther- 
malization at large distances. In other words, the Hawking radiation in AdS produces 
the dethermalization in dual QGP in the infrared region. To describe the Hawking 
radiation we use the AdS- Vaidya metric with negative energy (with negative mass 
function), compare with [f£3}| , where the Hawking radiation has been studied for the 
Minkowski Vaidya metric and see [B^] where the AdS- Vaidya with negative mass 
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has been considered. The process of the black hole Hawking radiation in this model 
corresponds to 

m(v) = M - M8(v), (3.1) 

and is presented in Fig.|l|. 



r-o 1= 



rt=rr 



X A. 



;:;:;:;:;:;;o;: : :: 



AdS 

-1/2 1/2 *B. 



:'jj;o:.jv(=i 




































.;j-BHAdSj_n 








i 



x c. 



r-o t>t Ihem , 

.-rrri-i-rrrfiSf.-rrrri--- 



X D. 



Figure 1: Cartoon of the black hole evaporation: A. The geodesic connected points ±1/2 
is totally in the black hole region. B. A partial evaporation of the black hole from the point 
of view of two points ±1/2 correlator, i.e. the geodesic is partially in the black hole region. 
. The total evaporation of the black hole from the point of view of this correlator, i.e. the 
geodesic totally abandons the black hole region. D. The geodesic is totally in the empty 
region in all subsequent momenta of time. 



The evolution of the 2-point correlation function ( |2.5| ) corresponding to the black 
hole evaporation process can be written explicitly in the case d = 2, see (|3~2"|), (|3.3| ) 
below. These formulas are similar to the evolution formulas describing the black hole 
formation [ 26| . The formulas are obtained by combining relations between the change 



of the value of the affine parameter and the value of the change of the coordinate x 
at the different parts of the geodesic. These relations have the form 



I — ^AdSA + lAdS,2 + iBHAdSA + ^BHAdS,2 



rl -Pi - 2p 2 x In ! " "r^A (3-2 



Px r c \ c *V \r*(2-r 2 G + + 2^F(nj) 

and 

Lren = 8L ren + LacLS,1 + LAdS,2 + ^BHAdS.2 + ^BHAdS,2 

Jl , 1 _ J?2 \ , r. 2 



-4 log(r c + Jr 2 r - pi) + In — v ; v 1 

V 1 -(^ + l-^) + 2r| + 2^R 



(3.3) 



Here ren means the renormalized length (compare with the action renormalization 
considered in |39|]) and the following notations are used 

G + = p 2 x + l-E B , G- = -pl + l + E 2 B (3.4) 
D(r) = r 4 + (E 2 B -l-piy+p 2 x (3.5) 
F(r) = P y -(pl + l- E\plY + 1 (3.6) 
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The above formulas correspond to the case when in the empty AdS space the geodesic 
has zero energy (the case of nonzero energy corresponds to non-equal time correlator), 
and the energy in the black space (in BHAdS space) is defined from the refraction 
condition 



E 



B 



pi- 



rn 



Here r c is the coordinate of the crossing point, p x is the angular momentum that is 
the integral of motion and has no a junction under crossing the geodesic (as opposed 



to the energy, that has according to (|3.7| ) a junction under a crossing of the shell), 
r t is the turning point of the geodesic. There is a relation 



t± 



(3.8) 



Let us clarify the meaning of ( |3.2|) , (|373|), (|3.7| ) and (|3.8| ). Here r c is a free 
parameter that specifies the position of the shell, p x is a parameter that does not 
change in shell moving, Eb is given by ( |3.7|) , and r t given by formula ( |3.8|) . Therefore, 
under r c and p x fixed, formulas (|3.2|) and ( |3.3|) give the relation between renormalized 
geodesic length L ren and the distance I. 

The process of the black hole creation and subsequent radiation within our model 
is presented in Fig. 0. The metric of this process is given by the formula ( j2.1| ) with 
m(v): 



m (y) = M6{v) - Mi6(v - v^, 



(3.9) 



The case of the total dethermalization corresponds to M — Mi and will be 
considered in what follows. If the distance i is much less then v\, then at the given 
moment of time the geodesic cross the shell no more than 2 times. This case is 
schematically presented in Fig.|| 



X A. 



























3- 










-E>!f 




















\ i 






















<V < 







• J_.J_.J_.J_.J_.J_.J_.J_.J_. 



X B. 



C. 



X D. 



Figure 2: Cartoon of the black formation and the subsequent black hole evaporation: 
transition from A to B shows the black hole creation; C shows the black hole evaporation 
and D shows the total black hole evaporation. 

As we can see from the figures, for < t < v\ the length of the geodesic 



connected the points (±//2, to, Zq) is given by formula from |p6|| . We assume that 



here t t h erm < v\. Then the length of the geodesic does not change in the interval 
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Uherm < t < V\, Fig.||B. Starting from v\ — t the geodesic crosses the shell with 
the negative mass (the magenta line in Fig.Q. C) at 2 points, and the length of the 
geodesic is given by formulas ( p.2[ ) and ( |3.3D . As it can be is seen from Fig.^. D, 
starting from t > v\ + t det the geodesic related the points (±//2, t , z ) occurs totally 
outside the black hole and the geodesic length will be given by formula 5C = 2 ln(£/2). 

For large £ the geodesic cross the shell at four points, see Fig|| and the length of 
the geodesic is given by a more complicated formula, but in any case at large times 
the length does not change and is given by SC = 21n(£/2), i.e. the thermalization 
does not occur. 
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Figure 3: Thermalization process does not occur for a long distance correlations: B. and 
C. show that the black hole does not have enough time to form, since the evaporation 
process has already started; A D as in Fig. 2 



4. Bounds on Thermalization and Dethermalization Times 

In this section we establish bounds on thermalization and dethermalization times in 
the AdS-Vadiya holographic model. 

Let us consider the AdS-Vadiya metric ( [2.1|) with 

m(v) = M(6(v) -6(v (4.1) 

where M > and v\ is larger than the thermalization time, v% > tther- 

Let us suppose that the vector I, characterizing the equal times two points on the 
boundary, at which the geodesic states and ends, has only one nonzero component, 
li = £, and denote J\ = J. The thermalization time for the correlation function at 
for these two points is given by 

f°° dr 

r therm -J j r2{1 _My ^ 

here J is the first component of the conserved angular momentum related with £ 

r°° dr 
£ = 2J . (4.3) 

J j r i / (r 2 _ J2) (1 _ m 
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From these formulas we get that 



I 

where m 2 = M/J d p d and F(m 2 , d) is given by 

roo dp 

2/; 



Tther ~-F(m 2 ,d), (4.4) 



F(m 2 , d) = £ — (4.5) 



Now, in the 4- dimensional Minkowski space, i.e. for ci = 4 for the function 
F(m 2 , 4) as a function of m, < m < 1, we have the bound 

0.39 < F(m 2 ,4) < 0.5, (4.6) 

see Fig.|j. Therefore we obtain the bound for the thermalization time (d=4) 

0.39 < < 0.5, (4.7) 



The similar bounds take place for other d > 2 [26, 55], ^9 . 

The dethermalization time, under assumption that V\ > T t h and for the same 
points on the boundary, is defined by the formulas (|4.2|) and ( [4.3| ) with M — 0. 
Since F(0,d) = 1/2 we have 

Tj f = I (4-8) 

Note, that ( [4.8[ ) does not depend on the space time dimension d. 

From ( [4.7[ ), ( [4.8| ) we obtain the following relation between thermalization and 
dethermalization times for observables at the same distance: 

0.78 < — < 1 (4.9) 

Tdet 

We see that this ratio is universal and does not depend on the distance between two 
points till the distance is less then v%. 

Therefore, the minimal ratio of thermalization time to determalization time, 
that can be realized in the d = 4 AdS-Vadiya model is 0.78. Increasing d one gets a 
possibility to decrease this ratio, see [|.B. As it has been noted in involving 
the nonzero chemical potential one increases the ratio of T t herm to £, and therefore, in 
our model, this increases Ttherm/ Tder- One can also try to add by hands an effective 
locking potential, for example, the quadratic one. This corresponds to a change 
(1 — m 2 1 p d ) — Y (1 — m 2 j p d + qp 2 ) in ( flBl) . This locking potential decreases the ratio 
F(m, g,4) = T therm /T der , see Fig.|.C. 

It is known that the experimental data on heavy ion collisions determine the 
thermalization time as Tther ~ 1 fm/c and the dethermalization time as T^t ~ 10 — 20 
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Figure 4: A. The plot of F(m 2 ,A) as function of m. B. The plot shows the dependence 
of the below bound of F(m 2 , d) on dimension d. d=2 corresponds to the red line, d=4,6,8 
to blue, green and magenta lines, respectively. C. The plot of F(m?,q,4) as function of q 
and vn? = 0.99. 



fm/c, so T t her/Tdet ~ 0.1 — 0.05 . It seems at the first sight that it is difficult to explain 
this data using our model. This is in fact so, if one thinks that thermalization and the 
dethermalization have to be happen at the same scale, but it is not so if the scales of 
thermalization and the dethermalization are different. For the thermalization time, 



the relevant length scale according |26| can be taken about £ ~ 0.6 fm, that is the 
thermal scale I ~ h/T for the temperature value T ~ 300 — 400MeV at heavy ion 
collider energies, and one obtains the estimate Ttherm ~ 0.3fm/c, which is smaller 
then the experimental data. 

One can fit better the experimental data using the scale I ~ 2 fm. One of 
possible explanations of this scale is a classical estimation of the distance between 
the nucleons inside the nucleus. One gets this estimation by taking into account that 
the radius of the nucleus of Pb is about rpt> ^7 fm, and in the sphere with this radius 
one can pack 208 (A=208 for Pb) balls with radius 



H r Pb « 1. 07 fm. (4.10) 

Here t\k is the Kepler number t]k = 7r/\/l8 ~ 0.74. From this consideration it seems 
natural to take I = 2r n ~ 2 fm as a typical scale of the thermalization. In this case 
one has T t h er m ~ 1 fm/c. 

For an estimation of the dethermalization time, one can take as the typical scale 
the size of the nucleus, i.e. Idet ~ ^fpb ~ 14 fm. Then one gets the dethermalization 
time Tdet ~ 7fm/c. Note that this estimation gets a low bound, since in the model 
we have the free parameter v\. For the ratio one gets 

T^L = ^ther hher- = q gg . 2 ^ Q Q ^ ^ 

T det 0.5 • ither <-det 14 

which is in agreement with the experimental data. 
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5. Conclusion 

In this paper we have used the holographic approach to study the processes of ther- 
malization and dethermalization in strongly coupled field theories relevant to the 
QGP formation in the heavy ion collisions. 

We show that the evaporation of black hole due to the Hawking radiation, which 
is modeled by the AdS-Vaidya metric with negative energy, leads to an interesting 
phenomena in the dual theory - thermalization is possible only at small distances 
and impossible in the infrared region. 

We have shown that for the ratio of the thermalization to the dethermalization 
time, considering at the the same length scale, has the universal bound and does not 
depend on the scale. The assumption that thermalization and dethermalization take 
place at the same scale does not fit well enough with experimental data. However 
it is more reasonable to consider thermalization and dethermalization at different 
length scales, since thermalization is local, while dethermalization is more suitable 
to relate with largest space scale of the system under consideration. It is obtained 
that the ratio of thermalization time to dethermalization time is equal to 0.05 that 
is consistent with experimental data. 

Acknowledgments. The present work is partially supported by the following 
grants: RFFI 11-01-00894-a (LA.), NSch-4612.2012.1. (LA.) and RFBR 11-01-00828- 
a (I.V), NSch-2928.2012.1 (I.V). 
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A. Hawking radiation and nonequal time correlators. 



Considerations in Sect. 3 concern to the case E = in the AdS space. Here we 
consider geodesic with nonzero energy in the AdS space. One of such geodesic is 
presented in Fig.[| This figure corresponds to the case when the shell has no time 
to reach the horizon. 



r=o 




Figure 5: A projection of the three dimensional picture in coordinates t, r,x of evolution 
of the shell describing the process for black hole formation. The picture contains a geodesic, 
related two point on the boundary with different times and with non-zero energy in the 
empty space. The green arrows —> show the time direction. The blue arrows indicate the 
time coordinates. The down thick line shows the shell position at the moment t AdS = t\, 
and the up thick line shows it at the moment t AdS = t%. A part of the geodesic with 
nonzero energy in the empty space is depicted by the dashed line and a part in the black 
hole region by the thick dashed line, geodesic abandon the plane t = const. The tick green 
segments show the time junctions under transition from AdS to BHAdS. 

In this case we have 



^AdS,! + lAdS,2 + I BHAdS, I + lBHAdS,2 



(pi - El) \ V r ? J (Pi- E l 

p , ln| (2- rfG+ + 2v ^)(2-,| G+ + 2^) | (A2) 
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and 



LacIS,i + L AdS: 2 + L B HAds,2 + L B HAds,2 + reg. (A. 4) 

= -2 In (n + - (pg - Ef^j (r 2 + y/r* - (jP x - E^ 

+ ^ In + r 2 + v^)^ (-± G+ + r 2 + v^)) 

-2 In (_±G + + T f + v /0(^ 

TAdS,l + Aii + T B HAdS,2 + T B HAdS,2 + At 2 + TA d S,2 (A. 5) 



- «i) \ V r? / 2 n + 1 r, 



+ ^ (, 1 + M + I + Il„^ 
(p£ - #1) V V r i / r 2 2 r 2 + 1 

£ , /lGLr 2 -G + \ £ , /lGLrf-G+ 

-arctanh , arctanh - 



2E B I 2 ^D{nj I 2E B 12 /D(^) 



E B \2 ^D{n) 



E B = E,-^-^rl + El- P l (A.6) 
The turning point is now in BHAdS and it is satisfied 



Here we use 



Aii = AAds^BHAdstin) ee ***** - ^ = -± In tl-il - I (A.9) 

2 ri + 1 ri 



At 2 = A flHAtS W(ri) = ^ S - = - + 1 In J^_Li (A.10) 

r 2 2 r 2 + 1 

It is clear that if ri = r 2 , then <5i = Aii + At 2 = 0. 
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